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WHEN DOES THE INTERGALACTIC MEDIUM BECOME 

ENRICHED? 

Benjamin D. Oppenheimer 1 , Romeel Dave 1 and Kristian Finlator 1 

Abstract. We use cosmological hydrodynamic simulations including 
galactic feedback based on observations of local starbursts to find a 
self-consistent evolutionary model capable of fitting the observations 
of the intergalactic metallicity history as traced by C IV between z = 
6.0 — > 1.5. Our main finding is that despite the relative invariance in 
the measurement of Q(C iv) as well as the column density and linewidth 
distributions over this range, continual feedback from star formation- 
driven winds are able to reproduce the observations, while an early 
enrichment scenario where a majority of the metals are injected into 
the IGM at z > 6 is disfavored. The constancy of the C IV observations 
results from a rising IGM metallicity content balanced by a declining 
Civ ionization fraction due to a 1) decreasing physical densities, 2) in- 
creasing ionization background strength, and 3) metals becoming more 
shock- heated at lower redshift. Our models predict that ~ 20 x more 
metals are injected into the IGM between z = 6 ^ 2 than at z > 6. We 
show that the median C IV absorber at z = 2 traces metals injected 1 
Gyr earlier indicating that the typical metals traced by C IV are neither 
from very early times nor from very recent feedback. 

1 Introduction 

Metals observed in the intergalactic medium (IGM) via quasar absorption line 
spectroscopy are thought to originate from some form of star formation-driven 
outflows as dynamical disruption is too inefficient (Aguirre et al. 2001). When 
and how metals are injected into the IGM are still a matter of debate. Observations 
of C IV in the IGM show surprisingly little evolution between z ~ 5 — > 2 (Songaila 
120011 Schaye et al. 120031 & Songaila 2005[) , despite this being the age of peak star 
formation. These observations motivate early enrichment scenarios by primeval 
galaxies and/or Populations-Ill stars at z > 6 where physical distances are smaller 
and the gravitational potential wells are shallower. 
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In Oppenheimer & Dave H2006f) (hereafter OD06), we explore a variety of 
galactic outflow models in Gadget-2 cosmological hydrodynamic simulations 
with the specific emphasis on modeling the C IV absorption observations between 
z — 6.0 — > 1.5. We are able to reproduce the main Civ observations including the 
relative invariance of f2(C iv) by implementing galactic-scale outflows based on ob- 
served local starbursts (Martin 2005, Rupke et al. I2005f) . Only a narrow range of 
wind parameters using the physically-motivated momentum-driven wind scenario 
(Murray, Quartet, & Thompson 2002J) are able to enrich the IGM without over- 
heating it at early times. Our models favor metal enrichment tied to concurrent 
star formation, as opposed to an early enrichment scenario with the majority of 
metals injected into the IGM before z — 6. Observations supporting this include 
the ubiquitous outflows seen from star-forming galaxies at high-redshift galaxies 
(Pcttini et al. 2001} and enhanced Civ and Ovi observed around Lyman-break 
galaxies at z ~ 2 - 3 (Adelberger et al. EESIEHOS- 

In these proceedings, we provide further arguments that the majority of metals 
in the IGM cannot be injected at very early times. Instead, we show the increasing 
IGM metallicity is balanced by an increasing C IV ionization correction resulting 
in observations of Q(Crv) that remain relatively invariant from z = 6 — > 2. We 
further explore the age of the metals in the IGM traced by C IV absorbers. 

2 Simulations 

We use the N-body+hydrodynamic code Gadget-2 fSpringel l2005|l with improve- 
ments including metal-line cooling and tunable superwind feedback as described 
in §2 of OD06. Briefly, the simulations we discuss here use an = 0.3 ACDM 
cosmology in a 16 h~ 1 Mpc box represented with 2 x 256 3 particles, giving a gas 
particle resolution of 3.9 x 10 6 M Q . §6 of OD06 shows that this sufficiently resolves 
most C iv absorbers. 

For simplicity, we focus on only two of the six outflow models of OD06, the 
"constant wind" (cw) and the "momentum-driven wind" (vzw) models. The cw 
model assumes w w i n d = 484 kins -1 and a mass loading factor r\ — 2, while the 
vzw model uses the relations i> W ind °c a and r\ oc cr _1 (tr being the galaxy velocity 
dispersion). 

Spectra were generated from the simulations using our software specexbin, 
which first calculates the physical properties along a line of sight and then applies 
ionization fractions from CLOUDY tables of an optically thin slab of gas with 
an impinging ionizing background field from Haardt & Madau ( 2001) including 
quasars and star-forming galaxies. The resulting spectra were analyzed using the 
AutoVP package to generate Voigt-profile fits from which line statistics could be 
derived and the summed fi(Civ) was calculated. 

3 Continual vs. Early Enrichment 

Figure n displays the evolution of fi(Civ) for a variety of observations and our 
cw and vzw models. A comparison to 17(C) indicates significant ionization cor- 
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Fig. 1. Evolution of f2(C iv) and fl(C), the mass density of C IV and metals, respectively. 
The constant wind (cw) model is overly efficient at enriching the IGM early, and creates 
an excessively hot IGM resulting in a much higher C IV ionization correction that increases 
too sharply at later times. The momentum-driven wind (vzw) injects fewer metals in 
the IGM, but does not require as extreme ionization correction due to less heating of the 
IGM from slower wind velocities at earlier times. The Ryan- Weber et al. data point is a 
lower limit. 



rections, usually in excess of 10 are necessary to derive the true metallicity from 
fi(Civ). The ionization correction increases toward lower redshift for three main 
reasons. First and most importantly for the vzw model, the physical densities 
decrease at lower redshift meaning that C iv traces metals at higher overdensities, 
which occupy less volume at lower redshift. Second, the ionization background 
strengthens between z = 6 — > 2 pushing the metals traced by C iv to slightly 
higher overdensity. Third, the metals become more shock-heated at later times 
due to metal-enriched winds colliding with increasingly faster infalling material 
toward low redshift. Metals are pushed into regions lower overdensity and higher 
temperature where carbon becomes ionized into C v and C VI. This third effect is 
greatest for the cw model causing its downturn in fi(Clv) at z < 3. 

New observations at z > 5.5 (Ryan- Weber et al. 2006. Simcoe 2006) may sug- 
gest even more C IV at high redshift than our models produce. The momentum- 
driven galactic superwind feedback prescription based on observations of local 
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Fig. 2. These three panels show the Civ absorption in a 100 kms~ slice of the vzw 
model (15 /i~ 1 Mpc across). The left and right panels show the metals at z — 6 and z = 2 
respectively. The middle panel represents the early enrichment scenario where the metals 
in our simulation at z — 2 are placed at z = 6, thus assuming metals passively evolve 
from z — 6 — + 2. The C IV ionization correction is much lower at z = 6 resulting in ~ 4x 
more absorption at z — 6 and an fi(C iv) that must decline by a factor of at least two 
from z = 6 — * 2, which is not observed. We can therefore rule out that the majority of 
metals were injected into the IGM before z — 6. 



starbursts and star formation rates may not be enough to enrich the high- 2 IGM. 
We do not rule out the possibility of another form of metal enrichment from the 
earliest galaxies and stars, although we cannot conclude this definitely. Neverthe- 
less, our point is to show that metal enrichment from continual star formation is 
required to produce most of the observed C iv by z = 2. 

As an exercise to show that the majority of metals cannot be injected into the 
IGM at early times, we take the metals in our simulation at z = 2 and put them at 
z — 6 densities and ionization conditions. This assumes that the metals passively 
evolve since z = 6 without any new metals being injected. We also allow the 
IGM to cool to the z = 6 IGM equation of state (EoS) at 5 < 100, thus assuming 
that the metals were injected early enough to allow the IGM to cool. Slices of 
absorption in Figure [21 show how this case (middle) would result in much more 
absorption than the original z = 6 or z = 2 vzw case (left & right respectively). 
The lower ionization corrections at z = 6 result in 4x more absorption leading 
to a predicted drop in fi(Civ) of at least 2x. Even if we leave the metals at the 
z = 2 shocked temperatures, there is still a significant decline. The fact that a 
constant Sl(Civ) is observed is strong evidence that metals in the IGM must be 
continually injected during the age of peak star formation. 

Finally, we explore the age of the typical C iv absorber at z = 2. We track the 
histories of individual SPH particles in our hydrodynamical simulation to see the 
redshift of their last injection into the IGM. Figure 3 plots the redshift histogram 
of a) the metal production (proportional to the global star formation rate), b) the 
last injection into the IGM of those metals, and c) the same as (b) but for metals 
traced by Civ. Although ~40% of the gaseous metals are in galactic/halo gas, 
only 15% of C IV absorption in spectra comes from this gas leaving the other 85% 
arising from the IGM. The median age of a C iv absorber (ie. the last time it was 
in a wind) is 1 Gyr at z = 2. Civ traces metals injected at a variety of ages, 
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Fig. 3. The black histogram (a) represents the fraction of metals produced at that 
redshift (proportional to the global star formation rate) indicating that 95% of metals at 
z = 2 are generated at z < 6 in the vzw model. On the same scale, the red short-dashed 
histogram (b) is the redshift when the metals were last injected into the IGM via a wind, 
while the blue long-dashed histogram (c) is the same except for metals traced by C IV. 
~ 15% of Civ absorption at z = 2 arises from galactic/halo gas, which accounts for 
~ 40% of the gaseous metals. The median C IV absorber traces metals injected 1 Gyr 
ago at 2 = 2. Thus C IV traces neither metals from only recent star formation nor early 
enrichment, but instead from a variety of ages. All histograms are on the same scale to 
allow one to compare them. 

and does not exclusively trace metals from very recent star formation nor metals 
from early enrichment. The average particle is injected into the IGM via a wind 
more than twice by z = 2, indicating there is significant recycling of metals in and 
out of the IGM. Older metals lie in lower overdensities suggesting observations 
that indicate C iv resides at a variety of overdensities (ie. the pixel optical depth 
method) require that metals be injected into the IGM at a variety of redshifts. 

4 Conclusions 

Our cosmological hydrodynamic simulations using momentum-driven winds from 
star- forming galaxies fit the main Civ observables between z — 6.0 — ► 1.5 as de- 
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scribed in OD06. We further emphasize: 

• The f2(Civ) measurement remains relatively invariant after z = 6 despite a 
vast majority (>90%) of the IGM becoming enriched between z = 6 — > 2. The 
ionization correction of Civ decreases due to 1) decreasing physical densities, 2) 
increasing ionization background strength, and 3) metals becoming more shock- 
heated at lower redshift. 

• We disfavor an early enrichment-only scenario, because it would lead to a de- 
creasing f2(Civ) with time due to the increasing ClV ionization correction. The 
IGM is primarily enriched via continual outflows from star-forming galaxies, al- 
though a different form of metal enrichment from the first stars/galaxies may be 
required for the z > 5.5 Civ observations. 

• C IV traces metals injected into the IGM at a variety of ages, with the me- 
dian age of 1 Gyr at z = 2. Only 15% of Civ arises from galactic/halo gas, 
while the other 85% is from the IGM. The spread of ages of IGM metals may be 
necessary to enrich the variety of overdensities as indicated by observations. 
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